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Abstract. We introduce the generalized rumor spreading model and investigate 
some properties of this model on different complex social networks. Despite pervious 
rumor models that both the spreader-spreader {SS) and the spreader-stifler (SR) 
interactions have the same rate a, we define a^^"^ and o;'-^-' for SS and SR interactions, 
respectively. The effect of variation of a'^^ and a^"^^ on the final density of stiflers is 
investigated. Furthermore, the influence of the topological structure of the network in 
rumor spreading is studied by analyzing the behavior of several global parameters such 
as reliability and efficiency. Our results show that while networks with homogeneous 
connectivity patterns reach a higher reliability, scale-free topologies need a less time 
to reach a steady state with respect the rumor. 
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1. Introduction 

Network modeling is the recent interest of a wide interdisciplinary academic field which 
studies complex systems such as social, biological and physical systems [H |^. By using 
a networked representation, it is possible to compare, in the same framework, systems 
that are originally very different, so that the identification of some universal properties 
becomes much easier. Moreover, a network description of complex system allows to 
obtain related information by means of completely statistical coarse-grained analyses, 
without taking into account the detailed characterization of the system. So universality 
and simplicity are two fundamental principles that are interested in the study of the 
emergence of collective phenomena in systems with many interacting components. 

From a general point of view, complex networks are connected graphs with, at 
most, a single edge between nodes where nodes stand for individuals and an edge 
corresponding to the interaction between individuals [31 H]. The collective behavior 
of nodes is complex in the sense that it can not be directly predicted and characterized 
in terms of the behavior of each individual. The collective behavior is the responsible of 
interactions that occurs when pairs of components are connected with links. It is simple 
to find various systems in both nature and society that can be described in this manner. 
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The most studied class of network modelling is the communication networks such as the 
Internet [5] and the World Wide Web [6]. A second class is related to social networks 
such as sexual contact networks [7], friendship networks [8] and scientific collaboration 
networks [9]. The last large class is concerned to biological networks such as metabolic 
networks [10] and food webs [H]. 

Rumors have been a basic element of human interaction for as long as people have 
had questions about their social environment [12]. Rumors are known for spreading 
between people quickly and easily since they are easy to tell, but hard to prove. 
Sometime, a rumor contains harmful information, so it is impossible to ignore, and 
can has damaging and perhaps even deadly consequences. We know it is probably bad 
for us, and we know it can hurt those around us, but we often find it hard to resist 
becoming active participants in the rumor spreading process. In business settings, it 
can greatly impact financial markets [T3| fT^. 

Despite its obvious negative connotations, a rumor has the capacity to satisfy 
certain fundamental personal and social needs and can shape the public opinion in 
a country [12]. To a great extent, rumors help people make sense of what is going on 
around them. In this case, rumors spreading becomes a means by which people try 
to get the facts, to obtain enough information so that it reduces their psychological 
discomfort and relieves their fears. 

A rumor can be interpreted as an infection of the mind. Daley and Kendall (DK) 
have introduced the original model of rumor spreading [151 IB]. In the DK model a 
closed and homogeneously mixed population can be classified into three distinct classes. 
These classes are called ignorants, spreaders and stifiers. The ignorants, those who 
have not heard the rumor yet, so they are susceptible to become infected by rumor. 
The second class consists of the spreaders, those who have heard the rumor and are 
still interested to transmit it. Eventually, the stifiers, those who have heard the rumor 
but have lost interest in the rumor and have ceased to transmit it. When the pairwise 
contacts between spreader and others occur in the society, the rumor is propagated 
through the population. If a spreader meets an ignorant, the last one turns into a 
new spreader with probability A; otherwise, the spreader meets another spreader or 
stifier, so they conclude that the rumor is known and do not spread the rumor anymore, 
therefore, turning into stifiers with probability a. An important variant model of DK 
is the Maki- Thompson (MK) model [17]. In the MK model, the rumor is spread by 
directed contacts of the spreaders with others. Furthermore, in the contacts of the type 
spreader-spreader, only the initiating spreader becomes a stifier. Therefore, there is no 
double transition to the stifier class. In the past, the DK and the MK models were used 
extensively to study rumor spreading [TH [191 120] • 

In the above-mentioned models of rumor spreading, the authors have investigated 
the rumor spreading in the homogeneous networks that their degree distributions are 
very peaked around the average value, with bounded fluctuations [211 122]. While, in the 
last years, a huge amount of experimental data yielded undoubtful evidences that real 
networks present a strong degree heterogeneity, expressed by a broad degree distribution 
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[221 121] • Recently, the model that we call the standard rumor model has been studied in 
Ref. |25] where authors studied a new model of rumor spreading on complex networks 
which, in comparison with previous models, provides a more realistic description of this 
process. In standard rumor model unlike previous rumor models that stifling process is 
the only mechanism that results in cessation of rumor spreading, authors assumed two 
distinct mechanisms that cause cessation of a rumor, stifling and forgetting. In reality, 
cessation can occur also purely as a result of spreaders forgetting to tell the rumor, or 
their disinclination to spread the rumor anymore. They took forgetting mechanism into 
account by assuming that individuals may also cease spreading a rumor spontaneously 
(i.e., without any contact) with probability 5. Furthermore, in the standard rumor 
spreading model, each node has an infectivity equal to its degree, and connectivity is 
uniform across all links. The generalization of the standard rumor model considered 
in Ref. [25] has been studied in Ref. [26] by introducing an infectivity function that 
determines the number of simultaneous contacts that a given node (individual) may 
establish with its connected neighbors and a connectivity strength function for the 
direct link between two connected nodes. These lead to a degree-biased propagation of 
rumors. To read more about social networks, one can refer to [271 l28l [29| 130] . 

In the above-mentioned models of rumor spreading, it has been assumed that in 
both interactions of stifling process (i.e., spreader-spreader and spreader-stifler) the 
initiating spreader becomes a stifler with the same rate a. In this paper we leave this 
assumption and assign a distinct rate for each interaction. More precisely, we introduce 
the generalized model in which the encounter of the spreader-spreader (spreader-stifler) 
leads to the stifler-spreader (stifler-stifler) with rate a^^^ (a^^^). We study in detail the 
dynamics of a generalized rumor model on some complex networks through analytic and 
numerical studies, and investigate the impact of the interaction rules on the efficiency 
and reliability of the rumor process. The rest of the paper is organized as follows. 
In Section 2, we introduce the standard model of rumor spreading and shortly review 
epidemic dynamics of this model. In Section 3 we introduce the generalized rumor 
spreading model and analytically study the dynamics of this model on complex social 
networks in detail. The influence of the topological structure of the network in rumor 
spreading is studied by analyzing the behavior of several global parameters such as 
reliability, efficiency in Section 4. Finally, our conclusions are presented in the last 
Section. 

2. Standard rumor spreading model 
2.1. Definition of model 

The rumor model is defined as follows. Each of the individuals (the nodes in the network) 
can be classified in three distinct states with respect to the rumor as /, the ignorant or 
the individual has not heard the rumor yet, the spreader or the individual is aware of 
the rumor and is willing to transmit it, and i?, the stifler or the individual has heard the 
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rumor but has lost the interest in it, and does not transmit it anymore. Based on Maki 
and Thompson model [17], the directed contact between spreaders and the rest of the 
population is the main requirement for spreading the rumor. From mathematical point 
of view, these contacts only can occur along the links of an undirected graph G{N, E), 
where N and E denote the nodes and the edges of the graph, respectively. The model 
that we call the standard model has been studied in Ref. |25]. By following [25], the 
possible processes that can occur between the spreaders and the rest of the population 
are 

• spreading process: SI — )■ SS whenever a spreader meets an ignorant, the ignorant 
becomes a spreader at a rate A. 

• stifling processes: 

— SS — )■ RS when a spreader contacts another spreader, the initiating spreader 
becomes a stifler at a rate a. 

— SR — > RR when a spreader encounters a stifler, the spreader becomes a stifier 
at a rate a. 

• forgetting process: — > R there is a rate S for a spreader to forget spreading a 
rumor spontaneously (i.e., without any contact). 



2. 2. Dynamics of standard rumor model 

The individuals in social complex networks not only be in three different states but 
also belong to different connectivity (degree) classes k, therefore we denote Ik{t), 
Sk{t) and Rkit) for densities of the ignorant, spreader, and stifier nodes (individuals) 
with connectivity k at time t, respectively. These quantities satisfy the normalization 
condition I^it) + Sk(t) + Rkit) = 1 for all k classes. We shortly review some classical 
results of standard model, where Nekovee et al. described a formulation of this model on 
networks in terms of interacting Markov chains, and used this framework to derive, from 
first-principles, mean-field equations for the dynamics of rumor spreading on complex 
networks with arbitrary degree correlations as follows: 

dikit) 



-kxikit)j2siit)pm (1) 

kXIk{t)J2Siit)Pm-kaSk{t)J2iSiit)+Riit))Pm-5Sk{t)i2) 



dt 

dSkit) 



= kaSkit) Y.iSi{t) + Ri{t))Pm + SSkit) (3) 

where the conditional probability P{l\k) means that a randomly chosen link emanating 
from a node of degree k leads to a node of degree /. Moreover, we suppose that 
the degrees of nodes in the whole network are uncorrelated, i.e., P{l\k) = lp{l)/{k) 
where p{k) is the degree distribution and {k) is the average degree. They have used 
approximate analytical and exacted numerical solutions of these equations to examine 
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both the steady-state and the time-dependent behavior of the model on several models of 
social networks such as homogeneous networks, random graphs and uncorrelated scale- 
free (SF) networks. They have found that, as a function of the rumor spreading rate, 
their model shows a new critical behavior on networks with bounded degree fluctuations, 
such as random graphs, and that this behavior is absent in SF networks with unbounded 
fluctuations in node degree distribution. Furthermore, the initial spreading rate at which 
a rumor spreads is much higher in SF networks as compared to random graphs. 

In standard model the authors have mainly focused on critical threshold in several 
models of social networks but in the following section we introduce generalized rumor 
model and we concentrate on the final fraction of the population that heard the rumor, 
R, when the spreading rate is fixed, A = 1, and we vary the value of other parameters 
of our model. 



3. Generalized rumor model 

In the standard rumor model [25], authors have assumed that both stifling processes, 
the SS — y RS and the SR — > RR, have the same rate a. But in this paper, we leave 
this assumption and define a^^^ and a*-^-* for SS — > RS and SR — > RR interactions, 
respectively. We will show that this separation of rates leads to notable results. The 
other interactions and their rates of our model are the same as the standard model. 
Now, the mean-field rate equations can be rewritten as 
dlk{t) k\ 



dt (k) 



hit)Y^lPil)Siit) (4) 



^ = ^^mT.mi)Si{t) - ^s,it)j:ipmit) 

ka^'^'^ 

- -— ^/P(/)i?Kt) - SSk{t) (5) 

^ = ^s,{t)j:ipm{t) 

+ ^s,it) y: ipmiit) + 6s,it). (6) 
\'^) I 

Eq. (4) can be integrated exactly to yield: 

hit) = 4(0)e-W^W, (7) 

where /a;(0) is the initial density of ignorant nodes with connectivity k, and we have 
used the auxiliary function 

= Yl kp{k) f Su{t')dt' = f{kSk{t'))dt'. (8) 

k "'0 "'0 

In order to get a closed relation for finding the final fraction of the population that heard 
the rumor, i?, it is more useful to focus on the time evolution of Assuming an 
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homogeneous initial distribution of ignorant, i.e., /^(O) = Jq (without lose of generality, 
we can put Jq ~ 1). The spreading process starts with one element becoming informed 
of a rumour and terminates when no spreaders are left in the population, i.e., S'fc(oo) = 
thus according to normalization condition, at the end of the epidemic we have 

Rkioo) = 1 - 4(oo) = 1 - e-W^(-) (9) 

After rather lengthy calculations, similar to what has been done in Refs. |25[ 126] , one 
can find the expansion of 0(oo) as 

' foo) = -^r^^ ttS — T (10) 



^ w 



1 + aW52M/2 + a(2)5/M(l - 51) 



where J is a finite and positive integral that has the form / = ^^^^ e^^^~^'^(f){t')dt' . At 
the end of epidemic, the final fraction of the population that heard the rumor, R, is 
given by 

i? = ^p(A;)(l-e~W<^(°°)). (11) 

k 

Regardless of the network topology and configuration, for any form of p(k), above 
relation can be simplified by expanding the exponential for the first order in (f){oo), one 
obtains 

i?~A0(cx)). (12) 



4. Results and Discussions 



One of the most important practical aspects of any rumor mongering process is whether 
or not it reaches a high number of individuals that heard the rumor. This value is simply 
given by the density of stifiers, R, at the end of the epidemic and is called "reliability" 
of the rumor process. For obvious practical purposes, any algorithm or process that 
emulates an effective spreading of a rumor will try to find the conditions that under 
these the reliability reaches as much as possible value. Another important quantity 
is the efficiency of the process which is the ratio between the reliability and the load 
imposed to the network. Load means number of messages on average each node sending 
to its neighbors in order to propagate the rumor. For these purposes, one does not only 
want to have high reliability levels, but also the lowest possible cost in terms of network 
load. This is important in order to reduce the amount of processing power used by 
nodes participating in the spreading process. 

In order to characterize this trade-off between reliability and cost, we use time as a 
practical measure of efficiency. Similar to Ref . [31] , we call a rumor process less efficient 
than another if it needs more time to reach the same level of reliability. 

To illustrate the effect of separation of stifling process rate a (for both SS and 
SR interactions) into a^^^ and a*^^^ for SS and SR, respectively, we consider a standard 
scale-free (SF) and Erdos-Renyi (ER) network. The SF network has generated according 
to p{k) ~ k~^, the number of nodes is = 4000 and the average degree is (k) = 8. 
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The ER network is a homogenous network that has the size = 4000 with (k) = 8. 
Throughout the rest of the paper we set A = 1 without loss of generahty and vary the 
value of a*^^^ and a'^'^K Fig. 1 shows the time evolution of the density of spreaders for 
different values of the stifling process rates, a^^^ and a^'^\ when the forgetting process 
rate is 5 = 0.5. 

We define two different models according to variation of a*^^^ and a^^^ as following 

• Model 1: a'^^^ varies with condition {a'^^^ < and a^'^^ = l}, 

• Model 2: a^"^^ varies with condition {a^"^^ < a^^^ and a^^^=l}. 

We have performed large scale numerical simulations by applying two stated 
conditions on SF and ER networks. Figs .1 (Fig. 2) corresponding to model 1 (model 2) 
illustrates, as expected, that the number of individuals who spread the rumor increases 
as the stifling process rate cx^^^ (ct''^'') decreases. In the cases in which the 
fixed, i.e., {a''^\a^'^^) = {a^'^\a^^^), the maximum value of spreaders in Fig. 1 (a) (Fig. 
2 (a)), the case in which < a^'^\ is greater than the corresponding values in Fig.l 
(b) (Fig. 2 (b)), the case in which ct^^^ < Q;d\ although the hfetime of spreaders in 
latter is greater. On the other worlds, when the ct^^ +q;^^) is fixed and < a^'^\ more 
individuals participate in spreading the rumor. On the other hand, in model 1 the time 
it takes for S{t) to reach its final value, i.e., no spreaders are left in the population, very 
slightly varies with different amount of a'^^\ but clear differences arise between the time 
that spreaders die out for different amount of a^"^^ in model 2. 

Fig. 3 (a) (Fig. 3 (b)) shows the final densities of stiflers for SF (ER) network. 
It obvious that in model 2 on both networks, the lower a^^^ leads to higher rehabihty 
(blue-solid curves). On the other hand, the time it takes for R(t) to reach its asymptotic 
value slightly increases when a^^) decreases, the clear differences arise for the two lowest 
values of a^'^\ Imposing the condition of model 1 on both networks, the red-dashed 
curves, illustrates that the lower a^^^^ leads to higher reliability but unlike the previous 
case, the time it takes for R{t) to reach its asymptotic value slightly decreases when a^^^ 
decreases as the inset figures show. 

Generally, from Fig. 3, after the comparison of the cases in which {a^^\a^'^^) = 
(ct*^^), a^)), we can conclude that the model 1 (blue-solid curves) leads to more reliable 
rumor spreading model but under condition of model 2 (red-dashed curves) the society 
reaches the steady state with respect the rumor in less time. 

To make the comparison between the number of stiflers in the SF and ER networks 
at the end of epidemic, we have plotted the Fig. 4. As shows this figure, in the ER 
network the number of stiflers at the end of the process is definitely higher than SF 
network, so the SF network appears less reliable. It results that ER networks allow a 
larger reliability R to this epidemic process. This is not straightforward and one may 
think that the existence of hubs in SF networks helps propagate the rumor. However, 
a closer look at the spreading dynamics reveals us that the presence of hubs introduces 
conflicting effects in the dynamics. While hubs may in principle contact with a larger 
number of individuals, spreader-spreader and spreader-stifler interactions get favored on 
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Figure 1. Time evolution of the density of spreaders in SF networks of size iV = 4000 
with (fc) = 8 for (a) model 1 (b) model 2. In model 1 (model 2) from below a^^'' (q^^^'') 
goes from 1.0 to 0.2 at fixed increments of 0.2. 




Figure 2. Time evolution of the density of spreaders in ER networks of size N — 4000 
with (fc) = 8 for (a) model 1 (b) model 2. In model 1 (model 2) from below a^^-' (a^^-') 
goes from 1.0 to 0.2 at fixed increments of 0.2. 



the long run. More precisely, it is very likely that a hub in the spreader state turns into a 
stifler before contacts all its ignorant neighbors. Once a few hubs are turned into stiflers 
many of the neighboring individuals could be isolated and never get the rumor. In this 
sense, homogeneous networks allow for a more capillary propagation of the rumor, since 
all individuals contribute almost equally to the rumor spreading. On the other hand, 
the SF network reaches the steady state with respect the rumor in less time than the 
homogeneous network (at the same condition). In this sense, SF network has a better 
efficiency. 
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time time 

Figure 3. Time evolution of the density of stiflers for (a) SF network (b) ER network 
with the same size N = 4000. The red-dashed and blue-solid curves are the stifler 
density corresponding to model 1 and model 2, respectively. The lowest curve in (a) 
and (b) is a common density for both models when a^^-* = a^^^ = 1. From below, the 
variables go from 1.0 to 0.2 at fixed increments of 0.2. The insets show maximum value 
of stiflers (the endpoints of curves) at the end of the epidemic. 




time 



Figure 4. Maximum value of stiflers at the end of epidemic for SF network (red- 
dashed) and ER network (blue-solid). The circles and squares show the maximum 
value of stifler with respect the time for model 1 and model 2, respectively. From 
below, the variables go from 1.0 to 0.2 at fixed increments of 0.2. 
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5. Conclusion 

In this paper we introduced a generalized model of rumor spreading on complex 
social networks. Unlike previous rumor models, our model incorporates two distinct 
rates for stifling processes. We have defined a^^^ and a^^^ for SS — > RS and 
SR — > RR interactions, respectively. Our simulations showed that in the condition 
(a^^^Q;^^^) = (q;^^\ q;^-^)), when a'^^^ is smaller than q;^^\ the society reaches the steady 
state with respect the rumor in less time. On the other hand, when a^"^^ < a^^^ , the 
higher level of reliability is obtained. This result is valid for both homogeneous and 
heterogeneous (scale-free) networks. 

By analyzing the behavior of several global parameters such as reliability and 
efficiency, we studied the influence of the topological structure of the network in rumor 
spreading. Our results showed that while networks with homogeneous connectivity 
patterns reach a higher rehabihty, scale-free topologies need a less time to reach a 
steady state with respect the rumor. 
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